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Abstract 

We present results of an experiment performed in Lake Baikal at a depth of about 1 km. The 
photomultipliers of an underwater neutrino telescope under construction at this site have been 
illuminated by a distant laser. The experiment not only provided a useful cross-check of the 
time calibration of the detector, but also allowed to determine inherent optical parameters of the 
water in a way complementary to standard methods. In 1997, we have measured an absorption 
length of 22 m and an asymptotic attenuation length of 18 m. The effective scattering length was 
measured as 480 m. Using {cos 9) = 0.95(0.90) for the average scattering angle, this corresponds 
to a geometrical scattering length of 24 (48) m. 
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1 Introduction 



Propagation of light in optical media is governed by two basic phenomena: absorption and 
scattering. In the first case the photon is lost, in the second case it changes its direction. 
The inherent parameters generally chosen as a measure for these phenomena are the absorp- 
tion length Xabs, the scattering length Xgd and the scattering function (or scattering tensor if 
polarization is taken into account) (3 {6). 

The determination of optical parameters for the water of deep lakes or oceans meets considerable 
problems. Typical values for Xabs in the window of maximum transparency are about 20 m for 
clearest water in deep lakes ||T], 0, ^ and 50 m for deep oceans ^, ^. Another important 
property of natural water is the strongly forward peaked scattering, with (cos 6) = 0.85- 
0.96 in oceans ^ or lakes [|], H, 3 and scattering lengths of 10 -50 m. Because of the long 
absorption and scattering lengths and the steep scattering function, devices to measure these 
optical parameters in clear water have to be extremely precise and well calibrated. 

In vitro measurements suffer from the comparatively short base length of the laboratory de- 
vices and from the fact that the samples of the medium can deteriorate before being analyzed 
in laboratory. Furthermore, for in vitro measurements it is rather difficult to perform long-term 
experiments and to observe seasonal variations of the optical parameters. In situ devices are 
difficult to handle, particularly for greater depths. On the other hand, only in situ measure- 
ments allow long term observations of the medium, including the implications of temporary 
changing mechanical, chemical, biological and thermal conditions. 

In the course of preparing a deep-underwater neutrino telescope in Lake Baikal, in situ op- 
tical measurements have been carried out over several years. We have constructed a special 
device (see |l^) and developed proper methods P, |, |^ to determine and to monitor the 
hydro-optical parameters of the water body at the site of Baikal Neutrino Telescope. The 
data obtained from these measurements do not only yield necessary input parameters for the 
operation of the telescope, but are interesting as well for environmental science. 

In this paper we present the results of optical measurements obtained in a way complementary to 
standard limnological (oceanological) methods, using a laser as light source and an underwater 
neutrino telescope as light detector. This kind of telescope represents a new class of giant 
instruments coming into operation at various locations in deep oceans and lakes and even 



in glaciers |T^. Being permanently operational over several years, underwater telescopes can 
perform, as a by-product and for calibration purposes, measurements of optical parameters as 
well as of water currents, bioluminescence etc. 

We have measured the spatial and temporal distribution of monochromatic light from a pulsed 
isotropic, point-like source. From the analysis of the spatial distribution, we derive the values 
for absorption length Xabs and for the so-called asymptotic attenuation length Xasm- From the 
time delay of the short light flashes arriving at the photomultipliers of the telescope, we derive 
Xeff, which is an effective parameter describing scattering. 

We describe the Baikal Neutrino Telescope in section 2, and the laser experiments to determine 
the optical parameters of the water in section 3. In section 4 we derive the expressions for 
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the light field as a function of the distance to an isotropic source. The experimental data are 
analyzed in section 5, and values for parameters characterizing absorption and scattering are 
derived. Section 6 contains the conclusion. 



2 The Baikal Neutrino Telescope 

The Baikal Neutrino Telescope |]13|, |14| exploits the deep water of the Siberian lake as a de- 



tection medium for secondary charged particles, like muons or electrons, generated in neutrino 
interactions. Neutrinos are etheric particles characterized by their extremely rare interaction 
with all kind of matter. They are supposed to be generated in cosmic particle accelerators like 
the nuclei of Active Galaxies or binary star systems. Shielded by kilometers of water burden, 
underwater telescopes aim to detect the rare reactions of neutrinos from these exotic sources. 

A lattice of photomultiplier tubes (PMTs) spread over a large volume records the Cherenkov 
light emitted by the relativistic charged particles. The light arrival times at the locations of the 
PMTs are measured with an accuracy of a few nanoseconds. From the time pattern, the particle 
trajectory can be reconstructed with a resolution of 1-2 degrees. In addition to the arrival 
times, the amplitudes of the signals are recorded. Typical amplitudes for particles triggering 
the array are in the few-photoelectron range. 




Figure 1: Location of the Baikal experiment 
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The Neutrino Telescope NT-200 is being deployed in the southern part of Lake Baikal (fig.l). 
The distance to shore is 3.6 km. the depth of the lake is 1366 m at this location. Fig.2 sketches 
the instrumentation of the site. Strings anchored by weights at the bottom are held in vertical 
position by buoys at various depths. The deployment of the detector elements is carried out 
in late winter, when the lake is covered by a thick layer of ice. Three cables {1,2,3 in fig.2) 
connect the detector with the shore station. Each shore cable ends at the top of a string {4,5,6, 
respectively). String 7 carries the telescope. A special "hydrometric" string 8 is equipped 
with instruments to measure the optical parameters of the water as well as water currents, 
temperature, pressure and sound velocity. The spatial coordinates of the components of the 
telescope are monitored by an ultrasonic system consisting of transceivers 9-14, and receivers 
along the strings. The relative coordinates of the components are determined with an accuracy 
of about 20 cm. 




Figure 2: Overall view of the NT-200 complex. 1,2,3 - cables to shore, 4j5,6 - string stations for 
shore cables 7 - string with the telescope, 8 - hydrometric string, 9-14 - ultrasonic emitters. The 
detail bottom left shows two pairs of optical modules (OMs), together with the electronics module 
controlling the OMs. Shown are two pairs directed face to face. 

The telescope NT-200 consists of 8 sub-strings arranged at the center and the edges of an 
equilateral heptagon. The light detection elements are optical modules (OMs) containing a 
hybrid phototube Quasar-370 with a hemispherical photocathode of 370 mm diameter. The 
time resolution of the Quasar-370 is 3 nsec for a single-photoelectron signal and improves to 
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about 1 nsec for large amplitudes. The OMs are arranged pairwise, with the two phototubes 
being operated in coincidence (see fig.2). If an event triggers at least 3 pairs within 500 nsec, 
light arrival time and amplitude of each hit pair {"channel") are recorded. Each of the 8 strings 
of NT-200 will carry 24 OMs. 

The time calibration is performed with a help of a nitrogen laser positioned just above the 



array |16]. The 1-nsec light flashes of the laser are transmitted by optical fibers of equal length 



to each of the OM pairs. 



3 The Laser Experiments 



For the two experiments described in this paper, an additional laser device was operated at 
various locations with respect to the PMTs of the neutrino telescope. Short light pulses were 
emitted by the laser, traveled through the water, and were detected by the PMTs during normal 
operation, when the standard trigger on muons crossing the array (3-fold coincidence within 
500 ns) was set. 



The laser module [jT6[ contains a pulsed nitrogen laser (wavelength 337.1 nm) which pumps a 
dye laser with an emission maximum at 475 nm. This wavelength well matches the point of 
maximum transparency at 490 nm, where Baikal water has an absorption length of about 20 m. 



Power supply, 
Micro computer 

Nitrogen laser 

Dye laser 
Attenuator 

Pressure vessel 



Isotropizer 

Figure 3: Schematic view of the laser module 



The light pulses generated by the laser have a length of about 0.5 ns (FWHM), thus being 
shorter than the time resolution of the photomultipliers. An attenuation disk, moved by a 
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stepper motor, is used to attenuate the light beam in five steps from 100 % to 0.3 % of the laser 
intensity of about 10^^ photons per pulse. Both lasers, the attenuator, power supplies and a 



micro-computer are contained in a glass pressure vessel |T8[ of 155 mm diameter and 1 m length 
(see fig. 3). After having crossed the vessel wall, the laser beam is directed into a hollow sphere 
with a diffusely reflecting inner surface. Through small holes in the sphere the light is able to 
escape, with a nearly isotropic distribution. Due to the large distance to the detector, the laser 
module can be considered as an isotropic point-like light source. 

The laser can be triggered to generate cycles of light pulses. Each cycle consists of 5 series of 200 
equally intense pulses. The intensities of consecutive series differ by a factor of 3 to 4. Laser- 
induced events fulflll the trigger conditions for muons crossing the array. They can be separated 
from those using the periodicity of the laser pulses (9.1Hz). The remaining background from 
muon events ranges from 10~^ for low light intensities to 10~^ for higher intensities. 

In the following, data taken with 2 conflgurations are analyzed: 
• 1995 experiment 

In April 1995, an intermediate version of the flnal array NT-200 was deployed. It carried 
72 optical modules mounted at four short and one long string (see flg.4). The laser 
experiment has been performed immediately after deployment of this detector {NT-72). 
The laser module was deployed from the ice surface and placed at several locations 20 
to 200 m away from the PMTs. Preliminary results from this experiment have been 
published in 



PMT pairs 
Laser 
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Figure 4: Schematic view of the 1995 configuration of the laser experiment. Each filled circle represents 
a pair of optical modules (channel). The laser is indicated by the open circle. 
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• 1997 experiment 

The 1997 laser experiment was performed after having deployed the first string of the 
1997 season. The string carried 24 PMTs, with the highest at 1059 m depth and the 
lowest at 1128 m depth. The laser was moved along the string at a horizontal distance of 
12 meters to the string with the PMTs, starting at a height close to the top PMT. The 
last of the 15 positions was 221 m below the first, at 1286m (see fig.5). 



In both experiments, the times and amplitudes of the fired channels were recorded if the stan- 
dard trigger on muons crossing the array (3-fold coincidence within 500 ns) was fulfilled. 



12.24 m-.- 
1065.6 m^ 



y- 



y- 



220.8 m 



*• 1059.39 m 



68.75 m 



1128.14 m 



y- 



y- 



1286.4 m^^ 



Figure 5: Schematic view of the the 1997 configuration of the laser experiment. Channels are shown 
expUcitly as pairs of OMs. The various positions of the laser are indicated by stars. The OMs face 
down, with the exception of channels 2 and 11 (upward arrows) which are directed upward. 
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4 The light field of an isotropic source in a medium 



The photon field $ (^R, Q, produced in a medium by a pulsed point-like source of monochro- 
matic light is a function of distance R from the light source, direction vector Q, and time t 
after emission. This function can be calculated by Monte Carlo methods, with Xabs, ^sct and 
f3 (6) as input parameters. In this section we derive analytical approximations for the photon 
field and compare it to the results of Monte Carlo simulations. 

In the special case of a medium without scattering one has 



$ [R, n,tj=F (R) 6 [nj 6{t- R/v) 
where F{R) is the photon fiux: 



F{R) = I $ (i?, n, t) dndt = e ^1 . (2) 

For a medium with scattering, F {R) can be written as 

where, in the general case, /i (R) is some intricate function of R. This function takes into account 
the increase of the photon path due to scattering. It was shown in ||20| that 1 < /i (i?) < 1.01 
for distances R < 0.6 Xgct and steep scattering functions ((cos 9) > 0.9). At larger distances, 
however, the contribution from scattering becomes significant. In the following we estimate the 
infiuence of scattering on the photon fiux at moderate distances from the source. 

Due to scattering, a photon does not travel along a straight fine but follows some polygonal 
path with random vectors {i = 1, 2, ...n) which are the trajectories of a photon between 
two successive scattering acts. Thus, 



(n) = Xsct, (4) 



in ■ fi+i) = Xl^t (cos 9) . (5) 

It turns out that under rather general assumptions about the shape of the scattering function, 
the following relation takes place: 

(cos 9i^i+k) = (cos 9)'' , (6) 
where 9i^ij^k is the angle between the vectors and rj+fc- In this case 
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R 



\ 



.ct^ n + 2 ^ (n — A;) (cos 9)^ . 

i=l 

The average length L of the photon path for n successive scattering acts is 



1=1 j=l 



(7) 



L = n ■ X 



set- 



For simphcity we consider the case of (cos 6) close to 1. With (1 — (cos 6)) being small, from 
from eqs. (|^ and (H) one gets 



+ (l-(cos 6)). 

For water with strongly forward peaked scattering we finally obtain 



(9) 



F{R) 



(10) 



We note that in this equation A^ct and (cos 9) are combined in an expression Xsct/{^ — (cos 6*)), 
resulting in a reduction of the number of independent parameters in eq. (|I^) from three to two. 
This agrees with the conclusions drawn in from the measurements of the photon flux of an 
isotropic source in artificial media as well as in natural water. For all these strongly forward 
scattering media it was found that F (R) can be described as a function of only two parameters, 
Xabs and P, where 



A 



abs 



X 



set 



Since for (cos 0) ^ 1 one has (0^) ~ 2 (1 — (cos 9)), eq.(p!OD can be rewritten in terms of the 
parameters Xabs and P (see also [|^, |2^ ) : 



F{R) 



AttR'^ 





e " 



(12) 



Our Monte Carlo simulations |2^ confirm the results of |21| which proved for an interval (9"^) = 
0.114 to 0.582 and for distances R < 5 Xabs that the error in F{R) resulting from the reduction 
of Xscatt and P{9) to P is smaller than than 1%. 

With the notation 



A 



eff 



A set 

- (cos 9) ' 



(13) 
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eq.(p!OD can be written as 



FiR) 



1+^ 



(14) 



The definition of Ae// according to eq.(|T3D coincides witli tlie one for tlie diffusion case provided 
(cos 6) is close to zero (i.e. quasi-isotropic scattering) and in the expansion in powers of (cos 6) 
only the first correction term is taken into account. We emphasize that in our case Ae// is just 



a convenient, artificial parameter (comparable, e.g., to the parameter P in eq. ([12D ) and should 
not be associated with an "effective" increase of the scattering length due to the anisotropy of 
the scattering. 



Although eq.(|lj) was derived here for strongly forward peaked scattering functions, similar 
considerations can be also performed for the case where (cos 6) is not close to unity as, for 



instance, in glacier ice which is exploited for the AMANDA neutrino telescope [12, 25 



Eq.(0) may also be considered as a particular solution of the transport equation for the small- 
angle approximation. In the framework of the small-angle approximation the general solution 
has the form (see e.g. p6 



47ri?2 sinh(i?VP/A, 



'^abs. (15) 



Eq.(|T^) can be formally found from the solution ([TSD by expansion with respect to the small 
parameter R^/P/Xabs- The opposite case, R^/P/Xabs ^ 1, gives 



FiR) 



2-nRXabs 



2lT 



(16) 



In fig. 6 we compare the results of our Monte Carlo simulations (histogram) with the predictions 
of eq. ([T^) (dashed line) and eq.([T5|) (solid line) for Xabs = 23 m and P = 0.1. The dotted line 
represents the result of eq.(0) with f^{R) = 1 (the case without scattering). One observes a 
perfect agreement of the approximation (plSj) with the Monte Carlo calculations over a wide 
region of distances. For not too large distances also eq.(0) agrees with the MC data. 

In next section we will apply the expressions obtained to the experimental data. 



10 



50 100 150 200 250 300 350 



R, m 



Figure 6: Intensity vs. distance as obtained from Monte Carlo calculations (histogram), eq.(12) 
(dashed line) and eq.([T5|) (solid line). In all three cases Aabs=23m and P = 0.1 have been used. The 
dotted line is the result of eq.(p|) with /i(i?) = 1 



5 Analysis of experimental data 

5.1 Absorption 

1995 experiment 

The absorption of light on its way to the photomultipliers can be determined from the measured 
amplitudes. For the 1995 laser experiment, 24 of the 72/2 = 36 channels have been included in 
the analysis. The combination of amplitudes from different photomultipliers requests a careful 
calibration. After amplitude calibration, and using the well-known angular acceptance of the 
photomultipliers, the mean amplitudes measured by each channel of NT-72 have been converted 
to light intensities. In fig. 7 the dependence of the photon density on the laser distance is shown. 
The points represent the measured signals from the optical modules facing toward the laser, 
for different laser positions and intensities. The data are well described by the relation 

F(i?)~-le-T (17) 
over a wide region of distances R. This is the approximation of equation (|l4l) for the case 
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R/ (SXeff) ^ 1. With Ae// ~ 480 m (see sect.5.2)[] and R covering the range 30 -150 m, the 
latter condition is satisfied. Therefore, A in eq. ([T7|) can be identified with the absorption length 
Xabs- The value of Xabs obtained from the experimental data of 1995 is (19.9 ± 0.2) m [l7]. 




20 40 60 80 100 120 140 160 180 200 

R/m 



Figure 7: Laser experiment 1995: Light intensity as a function of distance to the laser. The dots 
show the measured values, the solid line shows an exponential decrease with A = 19.9 m. If scattering 
is included, the dependence shown by the dashed line is expected. 



1997 experiment 

In fig. 8 we present the average amplitudes multiplied with R^ vs. distances R, for one of the 
channels of the first 1997 telescope string. The data points can be grouped along five straight 
lines, each corresponding to a given intensity of the laser. Contrary to the experiment of 1995, 
the analysis has been performed for each channel separately. This procedure avoids biases 
due to possible errors in the relative amplitude calibration of the PMTs. The large number 
of channels (i.e the multiple independent determination of Xabs) and the high precision of the 
measured amplitudes and distances results in small errors of Xabs- The absorption length is 
Xabs = 22 m, with variations of about 6%, dependent on the covered range in R. Most of this 
variation (about 5%) is due to the change of Xabs with depth. The errors of the measurement 
itself are estimated to be smaller than 1% (i.e. 0.2 m). We attribute the discrepancy between 
1995 and 1997 results to time variations of the water parameters. 

-•^ Actually, from the rather precise description of the data by an 1/R^ ■ exp{—R/X) dependence, a lower limit 
of about 300 m for Ae// can be estimated. 
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40 80 120 160 200 

R,m 

Figure 8: Laser experiment 1997: Logarithm of average amplitude (A) times squared distance 
versus distance R for one of the 12 channels of the first string. Curves are to guide the eye. Each 
curve corresponds to one of the five different laser intensities. Saturation of the curves at low R is due 
to the limited dynamical range, flattening at high R is due to the noise. Fits have been applied to the 
linear part of the curves. 



5.2 Scattering 

To investigate scattering we evaluated the delay in arrival times of the photons. According to 
eqs.(|^), (|T^), these times can be written as 

where is the velocity of light in water. Note that this approximation applies only for multiple 
scattering under small angles. Strongly scattered photons (e.g. due to Rayleigh scattering) 
arrive much later and would not coincide with one of the many forward scattered photons. This, 
however, is a condition for a channel to trigger, since its two OMs are operated in coincidence. 
Eq.(0) does not include effects of absorption. Consequently, in reahty the average time delay 
is smaller than given by expression (|l^): photons traveling a longer path are absorbed and 
contribute less to the average delay. Absorption is described by the exponential probability 
distribution of eq.(14). 

With Xabs determined from the intensity-vs. -distance curve, the effect of absorption can be 
subtracted from the experimental data as well as from Monte Carlo data. Fig. 9 compares the 
experimental delays to delay curves expected from different assumptions on Xeff- For both 
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experiment and prediction the effect of absorption is subtracted so that eq.(|T8[) can be apphed. 
The sohd hne gives the best fit and corresponds to Ag// = 480 m. The error of Ag// is estimated 
as 15%. 




50 70 90 110 

R,m 



Figure 9: Laser experiment 1997: Average arrival times {t) vs. distance R. Dots represent the 
measured values, lines give the expectations for different values of the effective scattering length X^jf- 
In both cases absorption effects have been subtracted (see text). 



The value of Ae// obtained in the present work can be compared with the resuhs of the direct 
measurements of Xsct and (3 {6) cited in [Q, ^]. There, with the help of the measured scattering 
function [3 (9), (cos (9)) ^ 0.95-0.96 was found. For A^ct values of 15-18m have been obtained. 
In 1997 (not pubhshed), we measured (cos(^)) ~ 0.9 and A^ct = 40-50 m. With eq.(p!3|), this 
resuhs in Ae// = 300 -450 m and A^// = 400 -500 m (1997), compatible with Ae// = 480 ± 70 m 
determined in the present analysis. We note that for (cos (9)) close to 1 the result is extremely 
sensitive to the exact value of (cos (6')). Also, scattering parameters in Lake Baikal are known 
to vary strongly with time (see |2^). This makes the agreement even more satisfactory. 

5.3 Asymptotic attenuation length 

As follows from eg. (p!6D , at large distances 

1 R 
F (R) ~ — ^aani ^ 

R 

where 

, _ Xabs 

1 + Y Aafes/ Ae// 
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(19) 
(20) 



is the so-called asymptotic attenuation length. Eq.(|T^) describes the asymptotic behavior 
of F{R). In the asymptotic regime, angular and spatial dependences of $ (^R, Q, tj can be 
factorized (see, for example, |]2B|) so that 

J $ [r, n, t) dt = F (R) p (n) . (21) 

Therefore, the measurement of Xasm is not sensitive to the orientation of the PMTs and can be 
performed independently for various directions. 

Although expression (IT^) follows from eq.([l5|) only for very large distances, our Monte Carlo 

R 

simulations show that the e ^'^"■^ /R behavior is already seen for R pa 100 m if Xabs ~ 20 m, 
Xeff < 500 m and the channel faces away from the light source. Photons observed by these 
channels are preferently very often scattered and form a light field close to the asymptotic 
propagation regime. Thus the value of Xasm can be evaluated from the data of the 1997 laser 
experiment where the maximal distance R was approximately 180 m. In fig. 10 the R dependence 
of the amplitudes for channel A^-2, looking away from the source, is shown. The experimental 
data behave exponentially with an exponent l/I/f ~ 1/(17.6 ± 0.5) m~^. 




100 110 120 130 140 150 

R,m 

Figure 10: Logarithm of the product of average amplitude {A) and R versus distance R in the 
"asymptotic" regime for channel 2, looking away from the laser. 

Comparing Monte Carlo simulations to eq. ([l5|) we found that in the general case 0.95 ■ < 
Xasm < 1-05 ■ and so we finally obtain Xasm = (17.6 ± 1.5) m. 

Since there is a certain overlap between the regions in which both eq. ([T9|) and eq.(|l^ describe 
the photon flux as a function of distance, one could try to estimate Ag// from Xasm and Xabs, 
using eq. (pO[) . Unfortunately, the final result is very sensitive to the exact values of asymptotic 
and absorption length (cf. what was said above on the estimation of Ae// from (cos (6'))). 
Therefore, the available accuracy of our experiments does not allow yet to determine Ag// in 
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this way. However, Xasm itself can be applied as a useful parameter for the description of the 
photon flux at large distances. 

6 Conclusions 

We have used a neutrino telescope triggered by laser light pulses to measure the optical prop- 
erties of clear, open water with high accuracy. For 1 km deep water in Lake Baikal, in 1997 
we determined an absorption length of 22.0 ± 1.2 m and an asymptotic attenuation length 
of 17.6 ± 1.5 m. Quite different to standard methods, the effective scattering length was ob- 
tained from nanosecond timing measurements. It was measured as 480 ± 70 m which under the 
assumption of {cos6) = 0.95 (0.9) translates to a geometrical scattering length of 24 (48) m. 

In future, we want to refine these measurements using a larger number of optical modules and 
measuring over an extended range of distances between optical modules and laser. Also, we 
plan a continuous monitoring of the optical parameters over the full year. 

This work was supported by the RFFI grant 97-05-96466. 
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